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Abstract 

We investigate chiral symmetry breaking and color symmetry breaking in QCD. The ef- 
fective potential of the corresponding scalar condensates is discussed in the presence of 
non-perturbative contributions from the semiclassical one-instanton sector. We concen- 
trate on a color singlet scalar background which can describe chiral condensation, as well 
as a color octet scalar background which can generate mass for the gluons. Whereas a 
non-vanishing singlet chiral field is favored by the instantons, we have found no indication 
for a preference of color octet backgrounds. 



1 Introduction 



Instantons, being pseudo-particles associated with tunneling processes, generate genuine 
non-perturbative effects in QCD. In the seminal work of 't Hooft P it was realized that 
they mediate an effective interaction between (light) quarks [HlllEllllElinilZllH]- This 
"instanton interaction" is attractive in the color singlet channel; hence, instantons pre- 
sumably play a role in the mechanism of chiral symmetry breaking O^J^J. In addition, 
they also provide for an interaction in color octet channels or in color triplet and sextet 
"diquark" channels. Mean-field computations based on a point-like instanton interaction 
have been employed as a central tool for investigations of color superconductivity at high 
baryon density f21ElElElE]; or for a description the baryon and meson spectrum and 
interactions in the vacuum in a Higgs picture with spontaneous color symmetry breaking 
PT| 13 ^] . This phenomenologically quite successful scenario requires a quark-anti-quark 
condensate in the color octet channel, giving rise to the question as to whether instantons 
support quark condensation in this channel. 

Symmetry breaking by a condensation phenomenon requires an interaction that lowers 
the free energy if condensates are formed. The bosonic condensates can be quark bilinears 
or even higher-order composites. In the case of instantons, a rich interaction structure is 
indeed provided: for A'^f light quarks, instantons typically induce an interaction between 
2iVf quark fields, which can be paired in many ways. This is one of the reasons why 
instanton-induced multi-fermion interactions have often been used as a starting point 
for investigations in the mean-field approximation. However, in the approximation of 
a point-like multi-fermion interaction, mean-field theory is ambiguous: by means of a 
Fierz transformation, the quarks can be grouped in different ways. For example, prod- 
ucts involving color non-singlet Lorentz scalars can be exchanged by products of color 
singlets in vector or tensor representations of the Lorentz group and vice versa. In view 
of this ambiguity, the relative strength between color octet and singlet channels remains 
undetermined, since the color octet channels can be completely removed or enhanced by 
suitable re-orderings fH]- Similar problems arise for the other colored channels used in 
the high-density computations. 

For further progress towards reliable computations, the Fierz ambiguity of the mean- 
field computation has to be resolved. This can be done in different ways. A first possibility 
explicitly includes the fluctuations of composite bosons after partial bosonization. Then, 
the dependence on the particular choice of bosonization (Fierz ambiguity) gets substan- 
tially reduced, as demonstrated by functional renormalization group techniques [20] • A 
second approach attempts to resolve the ambiguity by explicitly taking the momentum 
dependence of the instanton-induced vertex into account. Finally, we propose a third 
method in this article that avoids altogether the use of the multi-fermion vertex and 
rather computes directly the instanton contribution to the free energy in the presence of 
selected condensates. The various approaches have different strengths and shortcomings, 
and a reliable picture will probably only emerge by a combination of them. 

The advantage of a study of the momentum dependence of the instanton-induced 
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vertex is based on the observation that pole-hke structures which arise from the effective 
exchange of quark-anti-quark or quark-quark bound states can be associated to the given 
channel of the bound state. In contrast to a point-hke interaction, such pole structures can 
no longer be moved to another channel by Fierz reordering. Momentum-dependent vertex 
functions can be dealt with using functional methods, such as Dyson-Schwinger equations, 
A^PI effective actions, functional renormalization group (RG) or suitable combinations. In 
particular, we envisage the functional RG as a promising approach towards a quantitative 
study for the condensation phenomena at hand, for reviews see fTI\ |221 I2H1 12^- The 
computation of the flow equations involves only a narrow momentum range around a 
given renormalization scale k, thus reducing the impact of an incomplete knowledge of 
the detailed momentum dependence of the full propagators and vertices; see, e.g., [23 
EEl I2Z1 1211 12H] • Its application to the present problem requires an implementation for the 
non-perturbative sector of gauge theories, e.g. [211 EHl 12^1 IHOl IHH 1211 12I|, also employing 
bosonization techniques as developed in |321 IHBl I2II1 OH I21j, or NPI- and NPPI- flows 
as discussed in [231 ESI OHl I21| • In particular, the Fierz-type ambiguity of the present 
problem can be resolved within a 2PPI-effective-action approach, since all possible (local) 
fermionic pairings are effectively taken into account by this approach jl9j . 

In this work, we consider a more direct approach to instanton-induced color symme- 
try breaking by taking advantage of the following observation: possible condensates can 
be viewed as background fields that are coupled to quarks and gluons via Yukawa and 
gauge interactions. We concentrate here on scalar color singlet and octet condensates. 
In presence of a singlet condensate, all three light quarks become massive, thus influ- 
encing the weight of the fermion determinant in the instanton calculation. Additional 
octet condensates induce a mass split between an octet of fermions (here associated with 
the baryon octet) and a singlet. Furthermore, all gluons acquire mass through the Higgs 
mechanism. Both effects modify the instanton contribution to the free energy. In partic- 
ular, the effective condensate-dependent gluon mass acts as an effective infrared cutoff, 
strongly suppressing the contribution of instantons of size larger than the inverse gluon 
mass. Furthermore, the infrared cutoff stops the running of the gauge coupling such that 
the gauge coupling remains small for sufficiently large octet condensates, and perturbation 
theory becomes applicable. By computing the instanton contribution to the free energy 
in presence of the condensates, we get access to those parts of the effective potential that 
violate the axial U(1)a symmetry. Under the hypothesis that these parts dominate the 
octet dependence of the potential, we may try to draw conclusions if the minimum occurs 
for vanishing or non-vanishing octet condensate. Our computation of this response is 
based on two theoretical concepts: on the one hand, the full functional integral is evalu- 
ated in the semiclassical one-instanton approximation. On the other hand, the decoupling 
of massive modes is taken care of by a proper threshold behavior of the running coupling, 
as it is suggested by the functional RG. 

Our method needs assumptions how the non-vanishing singlet and octet condensates 
influence the masses of quarks and gluons. In practice, this is done by an ansatz for the 
effective action which describes the couplings of quarks and gluons to the color singlet and 
octet condensates. Apart from the restrictions imposed by color and flavor symmetry, the 
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details of this effective action are not known. This is one of the most severe restrictions on 
the quantitative rehabihty of our computation. Nevertheless, the qualitative features of 
mass generation for quarks and gluons can be captured in a simple picture. We consider 
here a local interaction with low powers of the condensates, in particular the chiral color 
singlet scalar aab with flavor indices a,b, . . . and the color octet scalar Xab,ij with non- 
trivial structure for flavor and color . . .). Our ansatz for the interactions between 
the condensate fields, quarks and gluons can be summarized in the following Euclidean 
effective Lagrangian |18j . 

£ = tZ^^^pi^t/jj + ^F.^^Fj,,^, (1.1) 

+Z^tr + ^.tr {SVta^a} 

1 + 7^ -j-l — 7^ 

-iZ^^lj^[{ha6ij + hxtj) — h {ha^Sij + hxji) — - — 

Here, we have included all power-counting relevant and marginal interaction operators as 
well as an effective potential for the background fields for completeness. In Eq. fll.lj) . we 
treat a and Xij as 3 x 3 matrices in flavor space and contract over the flavor indices of 
the quarks. A successful phenomenology of QCD based on an effective Lagrangian of this 
form has been worked out in ^3 ^1 IHHj • 

Obviously, the reliability of our conclusions will depend on whether the ansatz (jl.lll 
gives a qualitatively correct picture for the response to non-vanishing condensates. We 
therefore present a few additional arguments for its motivation. Associating the con- 
densates a, X with corresponding fermion composites ~ ipip, the interactions of the type 
specified in Eq. fll.l|) arise naturally from fundamental QCD, as can be studied with tech- 
niques developed in jS2l E31 1201 EHl Oil • particular, box diagrams of the type shown in 
Fig. n play an important role. In this work, we choose the viewpoint that these effective 
interactions are present in the dominant momentum region for the instanton contribution, 
being generated by U(l)A-preserving interactions, also partly at higher momentum scales. 
We do not attempt here to compute the parameters appearing in the effective action ()1.1|) 
except for the gauge coupling. For a qualitative study, we treat the Yukawa couplings 
h, h as well as the wave function renormalization factors Z^, Z^r, Z^ as free parameters. 

In presence of non-vanishing background fields a chiral symmetry is broken, whereas 
X acts like the Higgs scalar, giving masses to gluons and quarks. For the present purpose, 
it suffices to investigate in detail the following two directions in field space: 

]_ . 1 

These configurations correspond to the condensates of standard chiral symmetry breaking 
and a color-flavor locked combination of quarks and anti-quarks, respectively. In this 
background, all fermions aquire mass, 

M^ = ha + ^hx, Ms = ha - ^^X, (1-3) 



0-ab = 0-5ab, Xab,ij = ^Xi^iaSjb- -Sij5ab)- (1-2) 
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Figure 1: Box diagrams with fundamental QCD interactions (left two diagrams) generate 
effective (nonlocal) four fermion interactions (middle). Using rebosonization |22] these can 
be translated to (approximately local) Yukawa interactions interactions with propagating 
composite bosons (right). 

with a split between the octet mass Mg and the singlet mass Mi for x 7^ 0. The fermion 
determinant in the instanton contribution depends only on Mg and Mi. We use the 
freedom of scaling of the fields a and x to set /i = /i = 1. In this normalization, a 
and X are directly related to the masses. In the Higgs picture of the QCD vacuum, the 
expectation value for Mg should be associated with the mass of the lowest baryon octet and 
Ml with a baryon singlet, possibly A(1405), yielding [23 Mg = l.lSGeV, Mi = -1.4GeV 
or (Jo = 866MeV, Xo = — 2.08GeV. In our approach, we treat a and x ^.s free variables. 
The octet condensate in ()1.3p provides for an equal mass for all eight gluons. 



Here, g is the renormalized coupling taken at an appropriate scale. The x dependence of 
g will be discussed in detail below. Then remains the only undetermined parameter 
of our ansatz. The phenomenological ansatz of fHl ISH] associates Mg with the average 

1 /2 

mass of the lowest spin-one meson octet, Mg ~ 850MeV and suggests Z^' ~ 1/15. 

The paper is organized as follows. In Sect. |2l we discuss the various effects of the 
quark and gluon masses on the instanton integral. In the subsequent Sect. 01 we dis- 
cuss the asymptotic behavior of the instanton contribution to the free energy. In Sect, 
m we investigate which condensate backgrounds are preferred by the instantons. Our 
conclusions are presented in Sect. |S1 

2 Effective potential in one instanton approximation 

Consider a given background of scalar fields a and x, as introduced above. Our aim is to 
compute the instanton contribution to the effective potential for a and x in the presence 
of fluctuating quarks and gluons. For homogeneous a and X; the effective action F thus 
decomposes into 



where Q denotes the spacetime volume. The non-anomalous contribution f/o(o", x) con- 
serves the axial U(1)a symmetry and will not be computed here. The anomalous contri- 



Mg = Z'J^glx 



(1.4) 



T[a, x] = n U{a, x) = ^ {Uo{a, x) + f/i„st(fr, x) + U*,,,{a, x)) 



(2.1) 



4 



bution t/jnst is induced by configurations with non-trivial topology, mediating also U(1)a 
violation j^U]. We determine this part in semi-classical approximation based on instan- 
ton methods. In particular, we resort to the approximation of a gas of dilute instantons 
in which f/inst. can be expressed by an integral over the instanton size p and the prod- 
uct of gluonic (inch ghosts) and fermionic fluctuation determinants in a one-instanton 
background (see Appendix . 

1 

f/inst = dp exp(-87rV/(p)) Agi(x,p) det M^,,,-. (2.2) 

Here, the exponential factor reflects the classical action of the instanton, and Agi sum- 
marizes the contributions from gluons and ghosts in the instanton background. The last 
factor, with 

M^^ij = -pij + aSij + Xij, (2.3) 

represents the fermion determinant which is of central interest to our work. In particular, 
it contains the zero modes of the Dirac operator which are responsible for anomalous 
contributions and give rise to a strong a and x dependence even for small values of these 
fields. 

It is useful to decompose Eq. (j2.2|) into a factor Cz(c"?X)P) arising from the fermionic 
zero modes, and another non- zero-mode factor (^{a,XiP) that summarizes all remaining 
(classical, gluonic, fermionic) contributions: 

Uinst{(r,x) = - j c?pCn(o-,X,p)Cz(cr,X,p)- (2.4) 

All contributions have been studied frequently in the literature, beginning with the seminal 
work of 't Hooft P . As important new aspects, we include the color octet scalar and take 
the threshold behavior due to decoupling of massive modes into account. 



2.1 Lowest order in the background fields 

Assuming that a and x are small compared to all other scales, their main influence arises 
from the zero-mode contribution. In particular, the non-zero-mode factor does not 
depend on the scalar fields to lowest order. For an SU(A^c) gauge theory with A'^f flavors, 
Cn reads"^ [T], 

Here, Ds is a scheme- dependent constant. A discussion of the scheme dependence includ- 
ing the difference between massive and massless regularization schemes can be found in 
Appendix ^ Our scheme has been motivated by the functional RG which generically 

^In Appendix O we briefly review the contributions from the zero and the non-zero modes starting 
from results given in [Tll41|. Moreover, we use this appendix to introduce our regularization scheme. 
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provides for mass- dependent schemes that automatically account for a proper decoupling 
of massive modes. This is a convenient feature of our RG-inspired scheme; however, we 
observe no qualitative scheme dependencies of our results. For example, to zeroth or- 
der in the fields, Ds in our RG regularization scheme is given by (see |12l ESI HO] and 
Appendix ^ 

Dm = DMs = ^2^N^^^T)\{N, - 2)! e^P(-l-51137iVc + 0.29175iVf) = 6.005 x 10-^ (2.6) 

where the last equality holds for A^'c = A'f = 3. 

As discussed in the Appendices IHl O om RG scheme is constructed such that it 
matches the MS scheme in the small mass limit. It was demonstrated in jH] that the 
MS scheme gives satisfactory agreement with lattice data in the ultraviolet. Without a 
color- flavor mixing mass matrix (x = 0), the eigenmodes of p are also eigenmodes of 
and we are led to [H IH] 

Cz(p,cr,x) = ((detflavor(V^o(a,i)|M^,ij|^o(^, j)))su(3) =: -p'^'V{a,x), (2.7) 

where the inner angled brackets denote the scalar product of the zero modes ipo, and 
the outer angled brackets denote a group average over all possible directions for the 
instanton in color space. In the last step, we have separated off the simple p dependence 
~ p^f and defined the auxiliary potential V{a,x)- We have also used the persistence of 
(quasi-)zero modes in the presence of the regularization [B]. For x 7^ 0, ^ and do 
not commute in general, e.g. for condensates x with ()1.2j) . The eigenmodes of p and 
do not agree anymore for x 7^ 0- Therefore, strictly speaking, Eq. ()2.7|) does not hold in 
general. However, in leading order of an expansion in x and a it holds true, as shown in 
Appendix O Inserting these findings into Eq. ()2.4j) . we obtain 

^inst(a,x) = Via,x) [ dpp^'Xnip) =■■ CVia,x). (2.8) 



For A^f < 4, C is a finite number for physically admissible running couplings from the UV 
to the IR as discussed in Appendix lEl For small a and x? the potential V{a,x) carries 
all dependence on the scalar condensates. 

So far our discussion has made no use of a specific color or flavor structure for the 
background fields. Let us now specialize to the condensates specified in Eq. p.2|) . Using 
the gauge-group averages computed in [2] we find (see Appendix for details) 

^(,,,) . H- i.x= . ^ + ^.f (. - ^.). (2.9) 

In this crude approximation where U{a, x) = CV{<^, x) with ( being a field-independent 
constant, we observe two fiat directions, a = —-^^X and cr = -^X^ but no global min- 
imum. In fact, V{a,x) is unbounded from below, similar to the findings in jT^]. In the 
present case, this simply signals the breakdown of the approximation of small a and x- 
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Figure 2: Schematic plot of V{a, x) at fixed a > 0. Without higher-order corrections, 
the potential is unbounded from below (solid line) with a local minimum at x = 0. If 
the cutoff mechanism provided by the higher-order corrections is strong (dotted line) the 
global minimum remains at x = 0. However, if the suppression sets in only at rather large 
values of x (dashed line) we have a global minimum at x 7^ in addition to a local one 
at X = 0. 




Let us assume for a moment that the potential becomes stable beyond this approxi- 
mation or by the inclusion of Uo{(t, x) (cf- Eq- (EH)))- Then one might speculate that the 
first fiat direction, a = which is a line of local minima for cr < 0, characterizes a 

global minimum (the second fiat direction is not even a local minimum). However, in this 
case, the ratio = f = ~ 0.068 is far from the phenomenologically reasonable range 

r ~ 0.4 Since V is completely determined by the zero modes of the massless Dirac 
operator, this fiat direction will not be lifted by the inclusion of higher order corrections 
in the bosonic fields in the 1-instanton approximation, as long as the split into zero- and 
non-zero-mode parts remains justified. A similar fiat direction was also found in [TUj . 

Let us furthermore assume that, for instance, Uq induces a nonzero VEV for a. Since 
V{a, x) prefers a positive a, the resulting potential V{a, x) in the x direction looks like 
the solid line sketched in Fig. |21 The case of no color octet condensate, x = 0? then 
is a local minimum. For larger x? the higher-order corrections from the non-zero-mode 
contribution and the threshold effects will set in, stabilizing the potential in x direction. 
Now it is a dynamical question as to whether this stabilization sets in early, i.e., for rather 
small X; such that no other minimum is induced (dotted line). Or stabilization could only 
modify the region of large x (dashed) line, such that the x^ term of Eq. (|2.9p wins out in 
between and induce a color octet condensate. 

The second scenario of color octet condensate formation seems more difficult to be 
realized, since the ~ term and the ~ ax^ are of opposite sign and the coefficient of the 
X^ term is rather small. Unfortunately, the small coefficients in front of ax^ and x^ in the 
potential ()2.9p will limit even the qualitative reliability of our investigation. As an effect 
of the color averaging, the potential in the x direction is almost flat for a given value of 
a, in contrast to the pronounced potential in the a direction. For a given cr, the weak 
dependence of V on the "direction" x/c" could easily be overwhelmed by corrections in 
higher orders in a and x that are much more difficult to control. Despite this caveat, a 
quantitative analysis remains interesting and will be presented in the next sections. 
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Let us close this lowest-order consideration with the remark that a and Xi general, 
are complex fields. However, complex field values typically lead to large CP violation, 
making them phenomenologically unacceptable; this is the reason why we restricted our 
analysis to real field values. If a non-trivial phase between the octet and singlet conden- 
sates is favored in case of non- vanishing this may lead to an argument against the 
formation of color octet condensates in general. In order to demonstrate this point we 
assume for a moment that the effective potential for the relative phase between x ^i-nd a 
is dominated by the small field instanton contribution U ~ C,V{a,x)- Then, real positive 
values of a would be preferred due to the instanton contribution. This would in turn 
lead to a positive "mass term" ~ (cf. Fig. |5(b)| ), originating from the ay^ term in 
()2.9p : for imaginary x = this turns into —a\xy- Combining this with the small x^ 
term, the relative minimum of V^(cr, x) for fixed a > 0, |x| > would occur for a complex 
CP-violating x- Unfortunately, the impact of this observation is weakened by the very 
small coefficients of the ax^ and x^ terms arising in our approximation. The approximate 
fiatness in the x direction makes the potential infiuence of other effects large. In this 
context we observe that the [/(l)yi-conserving part Uq in ()2.1|) also contributes to the 
effective potential for the phase between x and a, for example with terms ~ cf*'^X^ + c.c. 
Only the common phase of x a-^id a is protected by the f/(l)yi-symmetry and is uniquely 
determined by the instanton part. 

2.2 Beyond small condensates 

As demonstrated in the preceding section, the instanton-induced effective potential can, 
in principle, support a mechanism for spontaneous color-octet condensation. Whether or 
not this mechanism is realized, however, requires a study that is valid for larger values of 
(T and X- The consequences of large condensates are twofold. First, the fermion masses 
are no longer small. This affects the non-zero-mode contribution Cn as well as the running 
of the gauge coupling. Also a mixing between zero modes and non-zero modes is induced. 
Second, a color non-singlet field gives an effective mass to the gauge fields, which again 
modifies the running of the gauge coupling (now the pure gauge contribution). In addition, 
it provides for an effective infrared cutoff for the p integration. 

2.2.1 Effects on the running gauge coupling 

Fermion and effective gauge boson masses exert an immediate influence on the running of 
the gauge coupling. For momenta smaller than the mass of a given quark or gluon degree 
of freedom, the corresponding fluctuations of this degree of freedom are suppressed. As 
a consequence, these fluctuations do no longer contribute to the running of the coupling. 
This decoupling of massive modes can directly be implemented in the (3 function for the 
running coupling, which we write as 
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where Mg and Mg are the gluon and the octet masses given in Eqs. p.3|) . p.4|) . and k 
denotes an RG momentum scale. The threshold functions /g,f(x) approach unity for small 
argument, /g,f(0) = 1, corresponding to the fact that the physical or effective masses play 
no role in the UV k oo. For large argument, i.e., for momentum scales k below a given 
mass, the threshold functions drop to zero rapidly, /g,f(x ^ 1) ^0, which implements the 
decoupling of massive modes from the renormalization flow. The threshold functions are 
not universal but regularization scheme dependent. For generic mass-dependent schemes, 
the threshold functions interpolate smoothly between the two limits.^ For the explicit 
computations, we set the threshold functions equal, lg{x) = /f(x) = l{x), and use 

This is a typical form for a threshold function, occurring in calculations based on the 
functional RG. Of course, the one-loop form used in Eq. ()2.10|) only serves as an example. 
A similar analysis of mass threshold behavior applies to any loop order and even fully 
non-perturbatively. We would like to stress that it is this threshold behavior where the 
additional free parameter enters via Mg, cf. Eq. fll.4j) . 

As a result of this decoupling mechanism, the effective running coupling is now field de- 
pendent, g{k, a, x)- Inserting this into Eqs. ()2.5j) . ()2.8|1 results in an additional field depen- 
dence of the effective potential. Qualitatively, the gauge boson mass weakens the increase 
of the gauge coupling. Owing to the exponential of the classical action ~ exp(— Svr^/gf^) 
in Eq. ()2.5|) . this leads to a total suppression of the instanton contribution. The fermion 
threshold behavior has the opposite effect due to the minus sign in the /3 function, reflect- 
ing their charge-screening nature. 



2.2.2 Effect on the instanton determinant 

The condensates give masses to fermions and gluons, hence the corresponding fluctuation 
determinants have to be evaluated for this massive case. Let us flrst consider the massive 
fermion determinant, i.e., the non- vanishing shift of the fermionic (non-zero) eigenmodes 
due to the effective fermion mass Eq. p.3|) . This problem has been solved recently using 
an efficient method to perform the mode sum HB] . The result interpolates smoothly 
between the analytically known small and large mass expansions |171I3H1- These calcula- 
tions have been performed with a color singlet quark mass m in the MS scheme. Here, we 
neglect the difference in the effect of the singlet and octet quark mass and approximate 

m = ^V|Mi|2 + 8|M8P. (2.12) 
o 

^For mass-independent schemes such as the MS scheme, threshold functions do not appear directly; 
but in order to describe the physics above and below a mass threshold adequately, theories with the 
correspondingly different particle content have to be matched at the mass threshold. This can equally be 
described by an effective threshold function which changes its slope discontinuously at a mass threshold. 
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Figure 3: The solid line gives the interpolating function K{x), smoothly connecting the 
small mass (dashed) and the large mass approximations EHl HE] • 



For our purposes, we have to adapt the results of jlHl UHl HE] to our massive RG regular- 
ization scheme, as derived in Appendix m and use the following interpolating function 



NiK{x) 



lndet'(-^ + m) 



RG 



-^iVf(i7(x) + ^) + In det \-p + m 



MS 
2 ^ ^4 



(2.13) 



ln(x) — ai + 



ln(a;) + ai — — asx' 



1 + a^x'^ + asx^ + a^x^ 

with X = pm. The function H{x) is defined in Eq. ()C.9|) . and 

ai = 0.792, as = 3.58, = 0.0842, = 0.00115, = 23.5, ag = 9.28, 

The primed determinant det ' in ()2.13p is that in the space of non-zero modes. As shown 
in Fig. |3[ this function interpolates smoothly between the small- and large-mass regimes. 



A similar behavior as for the fermion determinant is expected for the non-vanishing 
gluon mass in the gluon determinant. However, this effect is sub-leading, the dominant 
effect being the modification of the classical action at the minimum, see e.g. ^ IH]: for 
constant this gives a contribution to the classical action ASd = — 67r^Z^|xpp^ and 
therefore a factor of (cf. Eq. ()2.2j) ) 

exp(-87rV^7'(p)) \iix, p) ^ exp(-87rV^?'(p) - Gn'Zjxl'p') Agi(p) (2.14) 

in the integral ()2.8p . To summarize, the full inclusion of a and x in the fermion determi- 
nant and the Higgs-type of contribution to the classical action result in our final formula 
for the effective potential {Nc = Af = 3): 

Uinst = D^V{a,x) J dpp~'^(^-^^^ exp 1^-^^ -67r2Z^|x|V + 3A:(mp)^ , 

(2.15) 

where K is given in Eq. (j2.13j) . and D^q is defined in Eq. (j2.6p . Once the running of the 
gauge coupling is specified, e.g., using the one- loop form of Eq. ()2.10|) and identifying the 
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RG scale with the inverse instanton radius, k = 1/ p, we can investigate the landscape of 
the instanton-induced contribution to the effective potential for a and x- Fo^' fixed a, an 
additional x dependence arises from the explicit term ~ |xp in the "classical part", the 
dependence of g{p) on x and the threshold effect K{mp). Our approximation of K[mp) 
reflects probably only poorly the dependence on the ratio x/'^j W6 have also neglected 
the mixing between the fermionic zero modes and non-zero modes which would modify 
V{a,x). 



3 Asymptotic behavior of the effective potential 

In order to obtain a more analytic understanding of the effective potential t/mst, let us 
investigate its asymptotic behavior for the different regimes of small and large fields a 
and X- Of particular interest is the interplay between this asymptotic behavior and the 
running of the gauge coupling. As an important caveat, it should be kept in mind that 
our derivation of the effective potential is based on the semi-classical instanton gas ap- 
proximation. This approximation implicitly assumes that the one-instanton contribution 
is small, which translates into a small value of t/inst- Therefore, whenever a large asymp- 
totic behavior of t/inst is encountered, this may not necessarily reflect the true behavior 
but rather signal the breakdown of the instanton-gas approximation. 

Our derivation of Eq. (j2.15p so far made use of the speciflc one-loop running of the 
gauge coupling given in Eq. (j2.1(jp . Assuming that the functional dependence on the 
running coupling holds also in the general case, we use the form of Eq. ()2.15|) also for other 
theoretically or phenomenologically motivated running gauge couplings. For deflniteness, 
we will use gauge couplings with the following infrared [p oo) properties 



const 



p~*oo 



p' 



ln(p)|Pi°g 

iPpowcr 



, ^pcrt(p)0(^QCD 



flxed point 

logarithmic divergence 
power law divergence 
6(- — Aqcd)oo perturbative 



> ,(3.1) 



with positive constants p\og and Ppower] each infrared behavior will be adapted to show the 
same decoupling properties for massive modes as displayed in Eqs. fl2.1()|l and ()2.11|) . For 
simplicity, we assume that this IR behavior does not depend on the number of fermions. 
These running couplings and the corresponding resulting instanton densities at vanishing 
external flelds are shown in Fig. [7| in Appendix |El 

The effective potential for non-vanishing a and x is strongly influenced by fermion 
and gauge boson mass effects, respectively. Hence, the effective potential is expected 
to behave differently in the various directions of the a, x plane. For deflniteness, let us 
investigate four cases: (1) a small, x = 0; (2) cr ^ oo, x = 0; (3) cr = 0, x small and (4) 
o- = 0, X ^ oo- 
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The p integration can be performed analytically by splitting the integration domain 
into several intervals, each of them dominated by a different effect. For case (1), for in- 
stance, there is a UV regime (0, Ppert) where the running coupling is well approximated by 
one-loop perturbation theory without the influence of mass thresholds. In the consecutive 
interval, (ppert, f ) (with c being some constant 0(1)) the gauge coupling is dominated by 
non-perturbative dynamics, but fermion masses still do not play an important role. Fi- 
nally, in the interval (^, oo) the fermions are heavy compared to the scale 1/p, and pure 
gluodynamics dominates the running coupling. Neglecting logarithmic dependencies on 
a or X, we find. 



A^a^-Po + 52a'^'"t\ (2) : a^oo, x = ^ 

4JVcPpowcr /O ON 

Agx^f + B^x , (3) : = 0, X small ' > 



^4X'-*, (4): (7 = 0, x^oo 



where all terms with coefficients Ai arise from the perturbative interval p G (0,ppcrt), and 
those with coefficients Bi arise from the various non-perturbative intervals. Of course, 
only the Bi terms depend on the form of the running coupling specified in Eq. ()3.H) : 
in fact, Bi and S3 vanish for the "perturbative" gauge coupling of Eq. ()3.1|) . Further- 
more, Pq = yiVc — lA'^f is the one-loop coefficient of the P function with fermions, and 
p'^ = IMc denotes the one-loop coefficient for pure gluodynamics. Not surprisingly, the 
logarithmically divergent and the fixed-point coupling yield the same results on this level 
of accuracy; hence the parameter p\og of Eq. ()3.1|1 does not enter Eq. ()3.2j) . 

For A'^f > 4 + 2A'^cPpowor; the small-a and small-x behavior may significantly be modified 
compared to the perturbative expectation o"^', x^'- Here the possible dependence on the 
infrared behavior of the gauge coupling appears to be important. However, this simply 
reflects the fact that naive IR convergence in the p integration is lost for iVf > 4, as we have 
already noted before. In this case, the convergence is now restored by a combination of 
the suppression due to the finite fermion mass exp(— A'^fi^(mp)) and the suppression from 

the gauge coupling f at the expense of a direct dependence on the IR behavior 

of the gauge coupling. 

In the large-field regime, the potential grows faster than a^, x^ only if /5o < which 
corresponds to theories without asymptotic freedom. Conversely, for asymptotically free 
theories, the instanton-induced potential will not dominate over the non-anomalous part 
of the potential IJq which can be expected to exhibit a cr^, x^ growth for reasons of univer- 
sality. In view of the caveat mentioned in the beginning of this subsection, we interpret 
this result as a successful self-consistency check of the instanton-gas approximation. Fur- 
thermore, for A'^c = = 3 and fixed a, the instanton contribution vanishes f/inst ~ X~^ 
for large x- The instanton potential is therefore stabilized in the x direction. 
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Figure 4: Sections of the instanton potential Ui^st{o',0) (left panel) and ?7mst(0,x) (right 
panel) in units of Aqcd- The different large-field behavior is clearly visible: the potential 
in the positive a direction goes to — oo whereas it rapidly approaches for large 
Note the different scales of the potential itself for the a and the x direction. (Both plots 
use Z^'^^ = 1/15 and a gauge coupling approaching an IR fixed point = 100. The 
renormalization-scheme-dependent constant Ds is scaled out.) 



4 Numerical analysis of the effective potential 

Guided by the analytic knowledge obtained so far for the effective potential, let us study 
our full result for the instanton- induced effective potential t/inst(o', x) obtained numerically 
from Eq. ()2.15|1 . For definiteness, we use - as an example - a one-loop form for the gauge 
coupling modified such that it approaches an IR fixed point at g^^ = 100 in absence 
of condensates; this is in the ball park of IR results from RG flow equations [23 • Our 
conclusions remain similar for all other running couplings proposed in Eq. 1)3.111 . 

Let us first confirm the asymptotic behavior obtained analytically above along the a 
and X axes. As is visible in Fig. |4(a)| displaying U{a,0), the potential along the a axis 
is unbounded from below for a — > oo. In particular, for N{ = = 3, the resulting 
asymptotics of Eq. ()3.2|) yielding |?7(cr, 0)| ~ ctit is confirmed. Let us stress that the 
overall sign of the potential is negative for positive a, since the integral in Eq. ()2.15|) 
is always positive and the prefactor V{a, 0) = —a^ is negative. From this, we draw 
two conclusions: first, the instanton potential favors chiral symmetry breaking, but the 
value of the condensate is not determined by the instanton potential alone (at least in our 
simple one-instanton approximation). Second, the complete instanton potential f/inst(c", x) 
cannot have a global minimum since there is a direction in which the potential always 
decreases. 

Next we consider a pure x field. Figure |4(b)| shows that the instanton potential 
becomes flat rather rapidly for large X) as expected from Eq. ()3.2j) . Negative x are 
clearly preferred.^ From Fig. |4(a)[ we can read off the location of the minimum: Xmin ~ 

■^Let us stress that the relative sign of a and x is indeed important, since it changes the parity of 
some particles in the spectrum of the model Moreover, owing to the U(1)a anomaly, it is not clear 
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Figure 5: Effective potential f/inst(o"7X) for different fixed a. The left panel displays the 
cases of 0" = 0, 0.06, 0.128, 0.20Aqcd in green, black, red and blue, respectively. We 
observe that for values of a > (Tcrit ~ 0.128Aqcd (red curve) the trivial minimum in 
the X direction becomes the global minimum. The right panel shows the potential for 
realistic values of a = 1, 2,4Aqcd in red, green and blue, respectively (descending order), 
clearly demonstrating the absence of an instanton-induced color octet condensate. The 
parameters are chosen as in Fig. 0] and all numbers refer to units of Aqcd- 

4.26Aqcd ~ 1.4GeV for Aqcd = 330MeV. This results in an encouraging Mg ~ 350MeV. 

Nevertheless, it is important to observe that the potential in the x direction is rather 
shallow compared to the a direction, the relative height being ~ 10~^. This is a direct 
consequence of the relative prefactors in the potential V{a,x) of Eq- (|2-9P - The depen- 
dence of the absolute value of x on our remaining free parameter is strong, whereas 
it remains weak for Mg. This is in agreement with the expectation that, in the octet 
direction, the threshold effect of the gluon mass Mg is much stronger than that of the 

1 /2 

quark masses ~ x- For our quantitative results, we use the value =1/15 which is in 
the phenomenologically acceptable range [THl 

Finally, let us study the complete potential depending on both fields a and x- Even 
though we have already observed that the instanton contribution alone does not have 
a global minimum, it is nevertheless worthwhile to look for a local one. Such a local 
minimum indeed exists at (cr, x) ~ ("0.27, — 4.2)Aqcd with the absolute depth of the 
potential being U{a, x)\\oc. mm. ~ — 1.42_Ds ~ 0.009 in units of Aqcd (and using -D^jg ~ 
6 X 10~^). Since this minimum has a too small |cr| ~ 90MeV and is extremely shallow, 
it is not physically acceptable. Any generic non-anomalous contribution Uo is likely to 
remove this minimum. 

Since the local minimum is not acceptable and a purely instanton-induced global 
minimum does not exist, let us redo our analysis with one additional assertion: we assume 
that the non-anomalous contribution f/o(cr, x) to the effective potential supports scalar 
singlet condensation via spontaneous symmetry breaking for sufficiently strong gauge 
coupling (in fact, this has been shown to happen generically in QCD-like theories in [32]). 

whether the sign in the Yukawa coupUngs can be rotated away by a chiral transformation. 
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Figure 6: t^inst(f7", x) for negative fixed a. On the left panel, we clioose a = 
0, —0.06, — 0.12Aqcd (green, red and blue). The global minimum then is always at non- 
vanishing (negative) x- On ^he right panel, we depict the behavior for somewhat more 
negative a = —0.1, —0.27, —0.44, — 0.61Aqcd (green, red, blue and black). After increas- 
ing up to cr ^ — 0.27Aqcd (red curve) the minimum in the x direction becomes more and 
more shallow with decreasing a. The parameters are chosen as in Fig. |1] 

In order to introduce as few parameters as possible, we simply assume that Uq fixes a non- 
zero value of a, leaving the detailed form of Uq aside. Now, since Uo is of non-anomalous 
origin, its form reflects the full chiral symmetry (even though it can exhibit a symmetry- 
breaking minimum), implying its invariance under a —a. Hence, Uq does not prefer 
a particular sign of a. By contrast, f/jnst does prefer positive values of a as displayed in 
Fig.|4g 

The instanton potential for small positive a is depicted in Fig. |5(a)[ We observe that 
the global minimum for x is non-vanishing only for very small values of a. Beyond the 
critical value of cTcrit ~ 0.128Aqcd < 50MeV the global minimum is at x = 0. This holds, 
in particular, for realistic values of cr = 1 . . .4Aqcd, which is plotted in Fig. |5(b)[ For 
cr > Aqcd the dominant feature of the x dependence of f/inst is simply the vanishing of 
t/inst for large x which results in a relative minimum at x = 0. Figure p(b)| summarizes 
one of our main results, namely, that the instanton-induced potential appears incapable 
of giving rise to an octet condensate in the present instanton-gas approximation. 

Incidentally, the us study the potential also for negative fixed a. Even though positive 
values are clearly preferred by Ui^st, observation only constrains the modulus of a to be 
in the realistic range \a\ ~ 1 . . .4. Negative sigma indeed always give rise to a global 
minimum at x 7^ 0, i.e. supporting a color octet; see Fig El However, for realistic values 
of a, the instanton-induced potential is extremely shallow again, such that this minimum 
is likely to be washed out by the non-anomalous part (unless the latter is either fine-tuned 
or supports an octet condensate itself). 
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5 Conclusions 



We have calculated the one-instanton contribution to the effective potential in a back- 
ground of classical bosonic fields coupled to quarks and gluons. One field, a, has the 
structure of the typical singlet chiral condensate and the other, x, exhibits a color-flavor- 
locking structure, as conjectured for the scenario of spontaneous color symmetry breaking 
in the vacuum. Beyond leading order in the background fields, we have included effects 
of quark masses on the running gauge coupling and the quark determinant. In addition, 
the color octet condensate works as a Higgs field for the gluons, providing for an extra x- 
dependent contribution to the classical instanton action. We work in a massive regulator 
scheme which makes threshold behavior more transparent. 

For the realistic case of A'f = 3 light quark flavors, the instanton potential is unbounded 
from below for a pure singlet chiral condensate, favoring a non-trivial value of this con- 
densate and chiral symmetry breaking. Although there exists a local extremum with 
non-vanishing octet contribution, it is rather shallow and thus likely to be washed out by 
non-instanton effects. Moreover, this extremum has a nearly vanishing value for the chiral 
condensate, making it phenomenologically unacceptable. As the potential is unbounded 
from below along the a direction, a global minimum of the instanton potential alone is 
excluded. Other stabilizing effects are typically expected from the U(l)A-preserving sec- 
tor. Therefore, we have investigated if a color octet condensate is favored at fixed singlet 
chiral condensate. For realistic positive values of the chiral field a ~ IAqcd, X = is the 
global minimum in the color octet direction. For negative a with a similar absolute value, 
there exists a minimum with non-vanishing octet condensate. However, this minimum 
is also unnaturally shallow and thus will presumably be washed out by non-instanton 
effects. Moreover, the negative sign of a is disfavored by the instanton contribution. 

An analytic insight into the effective potential can be gained from its asymptotic 
behavior for large background fields. It is interesting to note that the behavior for very 
large fields is inherently connected to asymptotic freedom. For instance, large color octet 
condensates lead to a high-scale decoupling, such that any non-perturbative IR behavior is 
screened; in this large-x direction, the flattening of the potential can directly be related to 
the perturbative approach to asymptotic freedom. Another extreme example is provided 
by the case where asymptotic freedom is lost, e.g., owing to too many fermion species; 
in this case, the potential may go to — oo faster than the fourth power of the fields, 
preventing stabilization by a renormalizable non-instanton potential for the bosonic fields. 
On the other hand, we observe a strong qualitative dependence on the non-perturbative 
IR behavior of the coupling only for very small fields and more than four light fermion 
species: here, the increase of the instanton amplitude depends strongly on the infrared 
details of the gauge coupling and may be changed from the naive a^f, x^'- 

In conclusion, barring large higher-order, e.g., multi- instanton, effects, we find that 
chiral symmetry breaking is supported by instanton effects. On the other hand, the 
issue of color octet condensation remains inconclusive. In our approximation, we find no 
evidence that instantons favor a color condensation in the vacuum. 
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A The one-instanton approximation and the instan- 
ton gas 

Here we recall the computation of the instanton-induced effective action F within the 
dilute-gas approximation. In a dilute and weakly interacting instanton gas, the dominat- 
ing contribution to the generating functional is the one-instanton and one-anti-instanton 
contribution. It reads 

Zr = -^{Ui{a,x) + U;{a,x)) (A.l) 

where VL denotes the 4- Volume, and Ui is the effective potential corresponding to the 
product of functional determinants of the fluctuating fields in this background. The anti- 
instanton contributes with f/Ai(o", x) = Ui{a,x)- Within the dilute-gas approximation, 
the contribution of the |?T,|-instanton sector is given by Z,„ = ^ ^( . This leads to the the 
full amplitude 

n=0 n=0 

where we have normalized the zero-instanton amplitude to one. To lowest order in the 
bosonic background fields, this holds because we included the influence of the fields 
(masses for the fermions) only for the zero modes. However, in the absence of gauge 
fields, the Dirac operator has no zero modes. Beyond this approximation, this normaliza- 
tion corresponds to a modification of the non- instanton U(l)A-symmetric contribution to 
U{a, x)- From ()A.2|) . we read off the effective action in the classical background of a, Xi 

T = -\nZ = -Zi, (A.3) 

which serves as the starting point of our investigation in the main text. 

B Different regularization schemes 

In this appendix, we discuss how to switch between regularization schemes. We shall 
use a scheme which manifestly exhibits the decoupling of massive modes. It has been 
shown in [B] that topological effects persist within the RG framework used in the present 
work. This applies, in particular, to the existence of zero modes [B]. Moreover, at leading 
order one only has to take into account the explicit mass or regulator dependence. This 
amounts to using the well-known zero modes 
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For non-perturbative problems involving mass threshold effects, as they are induced 
by the background fields in our case, such schemes are highly advantageous. Following 
m 132], a change of the renormalization scheme can be understood by comparing two 
integrals. 



/i 



1 



2^2 ' 



dS 



(27r)4 (g2 + rn 



2^4 ' 



(B.i: 



in the different regularization schemes. Here, Ji appears in connection with the zero 
modes and I2 with non-zero modes. Therefore, we can easily keep track of these terms. 
For example, the Pauli Villars scheme gives 
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(21n(Apv) -21n(m)), (B.2) 
(21n(Apv) -21n(m)) =/f^. 



whereas we find in dimensional regularization 
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967r2' 



2 ln(/i) - 2 ln(m) - 7 + ln(47r) + 0(4 
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(B.3) 
(B.4) 



Comparing Eqs. ()B.2|) and ()B.4|) 
dimensional regularization read 

h : In(Apv) ^ 

J2 : In(Apv) ^ 



the substitutions for a change from Pauli- Villars to 
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Using Eq. ()B.5|) . it is easy to check that, starting from Eqs. ()2.5|) . ()2.7p . we can obtain 
the corresponding result in dimensional regularization as given, e.g., in ^J. 

Both schemes discussed so far are mass independent. This originates from the fact that 
for fixed cutoff /i or fixed dimensionality 4 — n, the integrals /^^ and I'^^^^ do not vanish 
in the limit m ^ 00. A mass-dependent regularization scheme should implement this 
decoupling: massive modes should not contribute to physics below the mass threshold. 
For m ^ fi the integrals Ii and I2 should become small and vanish in the above limit 
fixed, m ^ 00. This can be implemented by defining /f^*^ 
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If^ with 
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dkk-Hi"^ 
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2m2A2 + 3A^ 
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[B.6) 
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In fact, this is not an arbitrary definition, but receives motivation from various sources. 
First, it is very convenient to have /i = J2, since, when changing from the common Pauh- 
Villars scheme to our scheme it is not necessary to distinguish between the different 
contributions from Ji and l2- The main reason, however, is the simple form of the one- 
loop flow equation for the gauge coupling in Eq. ()2.10j) which results from this choice. 
Finally, a deeper reason for the choice is that it corresponds to a simplified version of a 
typical functional RG scheme regularization (for more details see, e.g. jSHl)- Indeed, 

I k J {2nY dk'''^'^ + + R,{q^)r ^ ' 

is a typical expression for Ji when one defines perturbation theory from a flow equa- 
tion with regulator function Rk{q^). In a consequent RG calculation several different 
threshold functions similar to l{uj) appear. For simplicity we put Ji = l2- This is not 
an approximation but simply an implicit definition of the related regulator function Rk. 
For computations beyond the present qualitative setting we suggest using an optimized 
regulator jSOl EI] 

Rk{q') = {k' - q') 9(1 - f^) , (B.8) 
and its upgrades suitable for momentum-dependent approximations 



From Eq. ()B.6|) . it is easy to find the relation between Pauli-Villars regularization and 
our scheme, 

In(Apv) ^ ln(H +It\^)=-2 ^^^^ + ^ ^ " I^W ^^"'^ 

We emphasize that Eq. ()B.9|1 depends on the cutoff A as well as on the mass m of the 
particle in question. For large m, the mass acts similar to the cutoff A. This implements 
the decoupling of heavy modes. 

Finally, we exploit the freedom of redefining the coupling constant at one-loop order, 
such that it absorbs part of the finite changes discussed above. 



^7r2 87r2 



Css- (B.IO) 



modifying the perturbative expression for only at order g'^. This is often used to simplify 
expressions, e.g., in the transition from MS to MS, or to ensure direct comparability 
between different schemes. We will use this freedom below in Appendix O to facilitate 
comparisons between our scheme and the MS scheme in which most results are given in 
the literature. 



C Assembling the instanton integral 

In this appendix, we put together all the various pieces of the instanton size p integral, 
taking care of our mass- dependent regularization scheme. At fixed instanton size and 
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using Pauli-Villars regularization, the following functions contribute to the renormalized 
integrand: 

87r2 

/ci = exp(-S'ci) = exp( ^—), (C.l) 

9 W 

/I A^f 1 \ 

/nonzero = exp ( " g^^c Hf^p) ' «(!) + y Hpp) + '2N{a{-) j , (C.2) 

4/471 \^^' 

= 7[m) 

/fcrmion ]y~ ; ( ^ • 

A'" 

4 7j-2{JVc-l) 

^ = ^(iVe-l)!(iVe-2)!' ^^-^^ 

where M^'' represents V^(o", x) as defined in Eq. p.7|l and reduces to m^' if the bosonic 
sources only lead to a simple mass term m. 

Here /d is the contribution from the classical action, /nonzero summarizes all effects 
from the non-zero modes, /gauge is the contribution from the gauge and /fermion from the 
fermion zero modes. N collects some normalization factors and the group averaging. 
Combining all these contributions, we find the well known result, 

iVn/. - ^PVP— (^) exp (-^ +/^olnM) (C.6) 



^Pvp-'+^' ( ^ ] exp 



in 



2 



where we have used the one- loop relation between g'^{p) and g^{p) in the last step, and 
D is given in the Pauli-Villars scheme as, 

Dpy = exp (^-a(l)-2(iVe-2)a(^)+2iVf«(^)^ = 1.1506, (C.7) 
= 2i? - 1 ln(2) - ^ = 0.1459, a(l) = = ^ ln(2) - ^ = 0.4433, 



R= Y^(ln(27r)+7) + ^5:^ = 0.2488. 

s=2 



We point out that we still have g'^{fi) in the prefactor of the exponential which is an 
artifact of the one loop calculation. This will be rectified by higher-loop orders where 
the 'bare' g^{fJ^) in the prefactor is replaced by its running counterpart evaluated at the 
scale p (also at one loop order less than the corresponding one in the exponential). Since 
replacing g^{p) g'^ip) is the main effect of higher loop orders (apart from possible 
changes in the factor D), we account for these prefactors as well as for the term in the 
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exponential by hand without further calculation. In this way, we already arrive at Eqs. 
(|2.5|) . (j2.8|) (at least naively in the Pauli-Villars scheme), 



(C.^ 



RG{M^p, m,p) = exp {^N,H{M^p) - '^N,H{mip)^ , H{x) = ^ Hl+x") 



The final task now is to change from the mass-independent Pauli-Villars scheme to our 
mass-dependent RG scheme. Starting from Eq. ()C.8|) . this is immediately done using 
()B.9|) . resulting in an additional multiplicative factor, 

3 + 2x2 

"4(1+X2)2' 

(C.9) 

where Mg and mf are the gauge-boson and fermion masses, respectively. 

Finally, we make a last change and define a modified RG scheme via Eq. ()B.10|) with 

c.„,^._,„(W^). (c.io) 

where Dy^ is defined in fj2.6|) . By construction, this establishes that our coupling con- 
stant is equal to the one-loop MS coupling, also including the constant in the instanton 
integral. But most importantly, we have not absorbed the mass-dependent contributions. 
Therefore, our scheme is still mass-dependent and provides for decoupling of heavy modes. 



D Calculation of the Zero-Mode Part 

The Dirac operator p in the background of an instanton has A'^f = 3 zero modes, being 
fiavor copies of a fundamental zero mode. We show that in leading order 

Cz(p,cr,x) = (detflavor(V^o(a,OI^V',iilV'o(^ j)))su(3) =: -p^'^(cr,x), (D.l) 

where 

{M^^ij)c,p = {{-pij)af3 + {crSij + Xij)Saf3) , (D-2) 

with mass matrix cr+x, where x introduces color-fiavor mixing. The eigenvalues of are 
A„((T, x) with eigenfunctions ijjn{c, x)- The color-fiavor mixing term x does not commute 
with p and the ipn are not eigenfunctions of p for x 7^ 0) leaving ()D.1|) a non-trivial 
identity. The mass matrix reads more explicitly 



Xab, 



(^Sab^ij + —7=X ( ^aiSbj - -Sab^ij 



Sa(3, (D.3) 



where we have already absorbed the Yukawa couplings into the fields and used the color- 
fiavor structure p.2|) for the condensates. For clarity we have explicitly written out the 
spin indices a, (3. Using the group averages in ()D.1|) |21E3, we arrive at 
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It is left to prove ()D.1|) . With trivial flavor structure, Xij = 0? the determinant factorizes 
trivially, 

det = Cz(p, (T, x) det 'M^ , (D.5) 

where det' stands for the determinant on the non-zero mode space. For non- vanishing 
X f!D.5|l holds up to terms x^^- This is shown in an expansion about the determinant of 
—Jp+a with eigenvalues A„(cr, 0) = A„(0, 0) + 0" and eigenfunctions ip{o', 0) = ipn{0, 0). We 
label the zero modes of pwith ipno, = 1,2,3 with eigenvalues A„(,(s = 0) = a. There 
is no term linear in x the only invariant is trflavorX = 0. The quadratic term has the 
structure (JU2{(j)x'^ with finite limit m(0). It is evaluated as 

ld!det{-p+a + sx)\^^, = lj2^d',\n)llXm+J2(9sXnm\^) J] ' (D-6) 



where the limit s = on the right-hand side of ()D.6|) is understood. The term proportional 
to d'^Xn has the coefficient Hm^n containing at least two of the eigenvalues A„o(s 



0) = (T of the zero modes ipno- Hence it only contributes to sub-leading terms like pcr^x^- 
The term proportional to (dgX)^ has, apart from sub-leading terms, one contribution 
proportional to H/^no mo -^^ removing two zero modes from the product. Thus we have 

M2(a) a^f-2x2 ^ ^N,^2 J2 (9,A„J(9,A„,J det'M^ + 0{a\') , (D.7) 

mo<no 

where det'M^ = H/-^^ the primed product involves only the non- zero eigenvalues of p. 
The s derivatives in (ID. 71) follow as 



dsXn{s = 0) = ds{iJn\{-p+ Cr + SxMn) = (V^n|xl^n) + A„(9^ ( | = {tpn\x\i^n) , (D.8) 

where we have used that {dsipn\{—p+ c")|'^n) = {ds4'n\4'n) Xn and {ipn\{—p+ cr)\ds4'n) = 
Xn{i/Jn\dsi^n) ■ We arrive at 

M2(0)X^ = (^nolxl^no)(V^molxlV'mo) det'M^ . (D.9) 

mo<no 

Equation ()D.9|) extends to general Mn(0)x" with n < N{. We are specifically interested in 
TVf = 3 with the remaining cubic term U3{a)a^^~^x^ = '^3(c")x'^, in leading order, 

^3(0)x' = (n^^^olxl^no) ) det'M^,. (D.IO) 

\ no / 

This proves (jET}. We can also directly use Eqs. ()D.9|l . (jD.10|l to compute the and x^ 
coefficients as group averages (m2)su(3)5 (w3)su(3) with the help of ^j. We arrive at 

\no=l I su(3) 



\no=l / su(3) 

leading to (^1^ . 
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Figure 7: Instanton density (right panel) for various types of infrared behavior for the 
gauge couphng (left panel). Red: g ~ const., green: g ~ ln(p), blue: g ~ p^, black: one 
loop behavior with = oo for ^ < Aqcd- Strong non-perturbative behavior is modeled to 
set in at roughly g'^ ^ 30. We note that the integral over the instanton size remains finite. 
More quantitatively, this holds if the plotted quantity on the right panel vanishes faster 
than with e > in the infrared. Most importantly, the constant C occurring in our 
lowest order approximation to the effective potential, Eq. ()2.8|) . remains a finite number. 

E IR running coupling effects in the lowest-order ap- 
proximation 

Here we demonstrate that the qualitative features of the instanton-induced effective po- 
tential to lowest order in the scalar condensates is largely independent of the behavior of 
the running coupling. This can be deduced from a study of the p integration in Eq. ()2.8p . 
For constant g^{p) = g"^ and A'f < 4, the p integral is infrared [p — *• oo) convergent, but 
has a (naive) UV divergence. The one-loop running removes this UV divergence, because 
of asymptotic freedom. Then, the integration kernel behaves as ~ p^^''"'"3^f~^(ln(p))^^'' 
for small p, rendering the integral convergent in this regime. In the infrared, the situation 
is less clear, since one-loop running is certainly not a valid approximation for the gauge 
coupling. Nevertheless, the restrictions on the behavior of (?^(p) for p — > oo are rather 
mild. Indeed, for positive and well-defined g^{p) there are no restrictions at all for A'^f < 4 
(massless flavors). We can even allow for a diverging coupling at a finite infrared scale 
Pdiv In this case, it is reasonable to assume that the coupling remains infinite for even 
larger distance scales, such that the integrand remains exactly zero for all p > pdiv In 
Fig. [7| we plot the integrand for running couplings with different infrared behavior. It 
is our main conclusion that all reasonable forms for the running coupling in the infrared 
imply a finite constant ( in Eq. (|2.8|1 . 
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